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We demonstrate the existence of robust bulk extended states in the disordered Kane-Mele model 
with vertical and horizontal Zeeman fields, in the presence of a large Rashba coupling. The phase 
diagrams are mapped out by using level statistics analysis and computations of the localization 
length and spin-Chern numbers C±. C± are protected by the finite energy and spin mobility gaps. 
The latter is shown to stay open for arbitrarily large vertical Zeeman fields, or for horizontal Zeeman 
fields below a critical strength or at moderate disorder. In such cases, a change of C± is necessarily 
accompanied by the closing of the mobility gap at the Fermi level. The numerical simulations reveal 
sharp changes in the quantized values of C± when crossing the regions of bulk extended states, 
indicating that the topological nature of the extended states is indeed linked to the spin-Chern 
numbers. For large horizontal Zeeman fields, the spin-gap closes at strong disorder prompting a 
change in the quantized spin-Chern numbers without a closing of the energy mobility gap. 



PACS numbers: 72.25.-b, 72.10.Fk, 73.20.Jc, 73.43.-f 

Introduction. Topological insulators (TI) are mate- 
rials characterized by robust properties against smooth 
deformations and disorder. Representative examples of 
TIs are the Chern insulators (CI) [l], quantum spin-Hall 
(QSH) insulators [2H5] and the strong topological insu- 
lators [6l [7] . Samples of these materials display robust 
conducting states at the edges while being insulators in 
the bulk. The edge modes are connected to robust metal- 
lic states residing in the bulk and away from the Fermi 
level [H [9] . The robustness of the edge modes and the 
bulk metallic states is a manifestation of a bulk topolog- 
ical invariant and, as such, there has been a sustained 
effort focusing on understanding the bulk states in disor- 
dered TIs [IMS. 

In two-dimensional bulk systems, the extended quan- 
tum states can survive localization only in extraordinary 
circumstances. In CIs, the bulk extended states carry a 
Chern number and, as such, they cannot be destroyed 
until these numbers are "annihilated" through collision 
with other extended states. In QSH insulators with time- 
reversal symmetry (TRS), it is believed that a similar sce- 
nario happens, but this time the bulk extended states are 
associated with the Z2 invariant [3 , which is protected 
by the insulating gap and TRS. Consequently, TRS is be- 
lieved to be crucial for delocalization in the bulk QSH in- 
sulators. The lack of a definitive answer is due to the lack 
of a complete theory of the Z2 invariant for aperiodic sys- 
tems (for recent progress see |2Ql[2T]). For example, the 
behavior of the Z2 invariant is not understood when the 
insulating gap is filled with dense localized spectrum, and 
the numerical approaches based on the twisted bound- 
ary conditions (BC) were not able to probe the diffusive 
regime, being limited to small systems [121 |T7]. Never- 
theless, analyses based on direct methods such as level 
statistics [19 and the transfer matrix approach [TTJ [15] 



have established beyond doubt that the bulk extended 
states are present in representative models of disordered 
QSH insulators with TRS. 

An alternative approach is provided by the spin-Chern 
invariants C± J0l|22]. While, for QSH insulators with 
TRS, C± and the Z2 invariants predict perfectly con- 
sistent phase diagrams [191 ESI [24] , they give conflicting 
predictions when TRS is broken. An early study [11] 
on the bulk extended states in disordered QSH insula- 
tors concluded that the protection against localization 
comes from both topology and TRS: Lack of any of the 
two will result in the immediate destruction of the bulk 
extended states. This conclusion was aligned with the 
predictions based on the Z2 invariant and was quickly 
accepted by the community. But recent studies on QSH- 
like phases with broken TRS revealed a more complex 
picture. Ref. [25 showed that robust gapless edge states 
can happen in QSH systems with TRS replaced by other 
symmetries. Yang et al. [26 studied the Kane-Mele 
model [2, 3 with a Rashba coupling and vertical Zee- 
man field, breaking both TRS and inversion symmetry, 
and found a topological QSH-like phase characterized by 
spin-Chern numbers that remain quantized until the bulk 
energy gap closes. The quantization of C± indicates that 
the topological order of the QSH systems is intact, and 
thus extended bulk states should exist in the systems, 
even when the TRS is broken. 

In this Letter, we establish the existence of such ex- 
tended states in disordered QSH models with strongly 
broken TRS. The phase diagrams of the models are ex- 
plored using the level statistics analysis [27 , the lo- 
calization length calculations ^28^ i29j and the spin- 
Chern invariants [TOl [22] . We show that the Kane-Mele 
model H] with a large vertical Zeeman field displays 
a nontrivial phase diagram in the {Ep^W) plane {Ep = 
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Covariance of level spacings 

FIG. 1. Covariance of the level spacings ensembles collected 
at various ^'s, VK's and for h = hzOz with hz = 0.1 and 0.2. 
The simulation was completed on a 40x40 lattice and with 
200 disorder configurations. The dotted line marks the 0.178 
covariance of GUE. The windows span horizontally between 
and 1. 



Fermi level and W = disorder strength) , with a topolog- 
ical phase characterized by quantized spin-Chern num- 
bers, surrounded by CI phases. Each phase is completely 
surrounded by lines of robust extended states. The situa- 
tion remains the same for horizontal Zeeman fields below 
a critical strength. Beyond the critical strength, the spin- 
gap closes at strong disorder, so that the extended states 
disappear without a closing of the energy mobility gap. 
Our results can be consistently explained based on the 
non- commutative theory of the spin-Chern numbers. 

The non- commutative spin-Chern Numbers. In a lat- 
tice model with many-orbitals per site |n,Qf,<j), where 
a {— ±1) and a represent the spin and other quantum 
numbers, let be the operator a^ln, a, cr)=cr|n, a, a). 
If P denotes the projector onto the occupied states of 
the Hamiltonian i7, one can consider the operator P&zP 
and its spectral projectors P± onto the positive/negative 
spectrum. In the presence of disorder, one can associate 
non-commutative spin-Chern numbers [30] to P± [22] : 

C± = 27ri E {tro{P± [ - P±], 2, P±]]}} , (1) 

where IE{ } = disorder average, tro = trace over the states 
at n=0 and x = position operator. C± remain quan- 
tized and invariant as long as = E {tro{P±ic^P±}} < 
00 [l9l|3T]. X± can be viewed as localization lengths and 
A± < 00 is enforced by a positive energy mobility gap 
and by the spin mobility gap in the spectrum of Pa^P. 
The localization length of P will be denoted by A^;^ . 

Eq. [1] can be numerically evaluated using only one 
(periodic-) BC and one disorder configuration {C± are 
self-averaging) [19]. This and the precise quantization 
conditions make C± very effective tools [18]. Note that 
in general ^A|;^ < A?_+A^, so A± < 00 implies A^;^ < 00, 
and A^;^ 00 implies A_ or A+ ^ 00. Therefore, when- 
ever C± are seen to take quantized values, one can safely 
conclude that the quantum states are localized at Ep. 
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FIG. 2. Normalized localization length Xe(M)/M as a func- 
tion of disorder strength W at various energies E, calculated 
for a long tube of 2M-site circumference, in the presence of 
vertical Zeeman fields. Symbols with different colors corre- 
spond to different sizes. The system length is taken to be 10^ 
sites. 



This allows one to map the regions of the phase diagrams 
that can harbor extended states, which are necessarily lo- 
cated in between the regions where C± assume different 
quantized values. 

Analysis of the spin- gap. Sharp predictions are pos- 
sible when the behavior of the spin-gap is well under- 
stood. For clean QSH models with and without TRS, 
it is known that the spin-gap stays open when the sys- 
tem crosses the QSH phase boundary where the energy 
gap closes [TO] |23l |26] . This remains true in the weak 
disorder regime and, as such, a sudden change of C± 
necessarily implies delocalization at Ep. For QSH mod- 
els with TRS, the spin-gap is known to stay open at all 
disorder strengths [19]. Now, for systems with and with- 
out TRS, from {Pcr^Pf = P [l - {i[crz,P]f) P we see 
that the spectrum of P&zP near is determined by the 
spectrum of i[az^P] near ±1, the edges of its spectrum. 
Using P = §c{C ~ with C encircling the energy 

spectrum below Ep^ one has 



P] = ^ §c{H - C)-i [a., H]{H - 0-HC . (2) 

Note that [cFz^H] is independent of (non-magnetic-) 
disorder, so using the standard localization estimate 



-s\rL — nn\/ Xe 



{s < 1, 



E{\{n\{H - E)-'\m)n 
= generic constant, A^; = the localization length at E)^ 
which can be explicitly proven for many models [32^ , one 
can establish: 



E{|(n|i[a„P]|m)|} <c,[a„if]^A|^e-^l"— , 

where O = the largest matrix element of O. Thus the 
averaged ^[a^, P] decays exponentially and remains small 
when [az^H] and A^;^ are both small (same for fluctua- 
tions). In such cases the spectrum of ^[(j^, P] is expected 
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FIG. 3. (a) Phase diagram of a QSH insulator with a verti- 
cal Zeeman field. The phase boundaries (solid lines) harbor 
extended states. The dotted line shows the phase diagram 
when Vr = 0. (b) and (c) show a computation of C± along 
the paths (1) and (2) (see panel a), respectively, completed 
on a 50x50 lattice and with 10 disorder configurations. 




Covariance of level spacings 
FIG. 4. Same as Fig. 1 but for horizontal Zeeman field. 



to be absent or localized at ±1. The conclusion is that, 
in general, whenever we see a sudden change in C± and 
[az^H] is small, Ep necessarily crosses an energy region 
where A^;^ is large or infinite. 

Heuristically, we can go even further and make the fol- 
lowing predictions: 1) The spin-gap is insignificantly af- 
fected by a vertical Zeeman term (in fact it is enhanced) , 
thus, in a QSH system with only vertical Zeeman field, 
a sudden change in C± signals the divergence of A^;^; 2) 
For a horizontal Zeeman field, there is a direct effect on 
il&z^ P] and the spin-gap is expected to decrease and then 
close as the field is increased. Nevertheless, the spin-gap 
is expected to remain open for fields below a threshold 
value or at moderate disorder, where a sudden change in 
C± will signal again the divergence of A^;^ . 

Numerical results. Our numerical simulations are 
based on the Kane-Mele model [2l |3] : 

H = + ^^SO E • {dkm X dnk)Cm 

{nm) {{nm)) 

Y1 4i^z • {O- X dnm)Cm + E 4(^n + h • &)Cn , 
(nm) n 

where last term represents our addition of a spin inde- 
pendent, on-site potential Wn with random amplitudes 
uniformly distributed in [— VF/2, VK/2], and of a uniform 
Zeeman field that breaks the TRS. We will fix Vso and 
Vr at 0.1, so that in the absence of the last term the 
model is inside the QSH phase. 

Let us consider first a vertical Zeeman field h = hz^z- 
In Fig.[l]we show the level statistics analysis for hz = 0.1 
and 0.2 and various disorder strengths. The ensembles 
of energy level spacings were recorded from small energy 
windows centered at various energies E and their covari- 
ance is plotted in Fig. [l] as a function of E. For both /i^'s, 
one can see energy regions where the covariance is large, 
indicating localization, but also two regions where the 
covariance is practically equal to 0.178, the covariance of 
the Gaussian unitary ensembles (GUE), indicating the 
presence of extended states. Furthermore, the extended 
states are seen to drift towards each other and ultimately 



collide, annihilate and disappear as the disorder strength 
is increased. This is a strong indication that the extended 
states carry a topological number. For a system in the 
unitary class, one would expect extended states to occur 
at sharp discrete energies. This is precisely what was 
observed for CIs ^3 |^, but here we see in Fig. [l] en- 
tire energy intervals where the covariance stays at 0.178, 
signaling a more complex phase diagram. 

The delocalized character of the states can be firmly 
established by using the recursive Green's function ap- 
proach [28^, 'W . In Fig. |2]we plot the calculated normal- 
ized localization length \e{M)/M as a function of W 
for long tubes with circumferences ranging from M = 6 
to 32, and for hz = 0.05, 0.1 and 0.2, and E=-2, -1 
and 0. For each E, we can clearly see instances where 
Xe{M)/M decays sharply with M, but also where it re- 
mains nearly independent of M, essentially forming a line 
of fixed points. In the later case, Xe{M) grows linearly 
with M therefore diverging in the thermodynamic limit. 
Consequently, the states are delocalized at these E^s and 
W's. Furthermore, the independence of Xe{M)/M on M 
indicates that we are dealing with critical states charac- 
terized by a vanishing P scaling function [28", ^29^ , a char- 
acteristic of the Kosterlitz-Thouless type transition [33 . 
The extended states must have a topological origin, oth- 
erwise the system will belong to the trivial unitary class 
where all electron states are localized at positive VK's. 
The data in Figs. [T] and [2] are consistent, for example, 
for £^ = —2 we see the covariance staying close to 0.178 
even at = 3 and after that the extended states drift to- 
wards E = 0. At about the same W, we see Xe=o{M)/M 
changing its behavior in Fig. |2|cl). The extended states 
are seen to collide at about W = 5 and £^ = in Fig. [ij 
and this is exactly where we see Xe{M) diverging in 
Fig.[2]:a3-c3). 

To construct the phase diagram, it is convenient to 
start from the limit Vr = hz = 0^ while holding Vso 
at 0.1. In this limit, the spin up/down sectors decouple 
and the phase digram in the {Ef^ W) plane, for each sec- 
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FIG. 5. Same as Fig. [2] but for horizontal Zeeman field. 

tor, consists of CI phases (with C± = ±1) surrounded 
by a line of extended states and then by a trivial insu- 
lating phase [ini [19]. TRS assures that the phase di- 
agrams for the spin up/down sectors overlap perfectly. 
As in Fig. [3|a), when hz is turned on, the phase dia- 
grams shift revealing two regions where the total Chern 
number C = C_ + C+ takes the values ±1, surround- 
ing a region which can be characterized by C± = ±1. 
The stability of the Chern numbers assures that the CI 
phases do not disappear when Vr is slowly turned on 
but they could part away and open the phase diagram 
at large W's. This cannot happen here because the spin- 
gap stays open so the C± = ±1 and C± = phases must 
be separated by extended states, a fact that can be seen 
explicitly in Fig. |2|a3-c3). We must conclude that the 
C± = ±1 phase remains completely surrounded by the 
CI phases as illustrated in Fig. [sj^a). The calculated val- 
ues of C± in Fig. [sj^b-c) for hz = 0.2, along the paths 1 
and 2 in Fig.jsj^a) further confirm our conclusions. Along 
path 1, both C± remain strictly quantized to ±1 until 
they simultaneously start an abrupt descent to 0. C± 
are seen to cross the values ±0.5 at about = 5, ex- 
actly where the annihilation is observed in Figs, [l] and |2j 
Along path 2, C- remains quantized at —1 well after C+ 
started its descent to 0, indicating indeed the crossing of 
a region with total Chern number -1 (also confirmed by 
a direct calculation of C). This also provides direct evi- 
dence that the spin-gap remains open when crossing the 
boundary between C± = ±1 and the CI phases, because 
its closing would have effected the quantization of both 
C± but in Fig.[3];c) C_ clearly remains quantized. A last 
remark here is that, since the CI phases are surrounded 
by C = phases, they are necessarily surrounded by a 
line of extended states as illustrated in Fig. [sj^a). 

We now consider the horizontal Zeeman field. In Fig.|4] 
we show the level statistics analysis for h = hx^x with 
hx = 0.1 and 0.2. At smaller W's we can again see 
regions where the covariance stays very close to 0.178, 
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FIG. 6. Same as Fig. [3] but for horizontal Zeeman field. The 
solid lines in panel (a) harbor extended states. The dotted 
line is where the spin-gap closes. 



and we still see these regions drifting towards each other. 
However, the covariance is clearly seen to move away 
from 0.178 before the annihilation is complete (see pan- 
els a5 and b5). This is a first indication that the line of 
extended states is broken, which is indeed confirmed by 
the recursive Green's function analysis shown in Fig. [5j 
At £^ = for hx = 0.1 and 0.2, in panels (b3) and (c3) 
where the annihilation more or less takes place, we see 
Xe{M)/M decreasing with M, indicating a finite local- 
ization length. Nevertheless, both Figs. |4] and [5] indicate 
that extended states are still present at lower W's. This 
is also evident in Fig. 2 of Ref. [11 . Furthermore, the 
data in Fig. |5] for the lower hx = 0.05 clearly demon- 
strate that the extended states survive here all the way 
till the annihilation point. 

Based on the numerical data, on the additional com- 
putations of C± in Fig. [6jb-c) and on what we learned 
from the case with vertical field, we derived the phase 
diagram illustrated in Fig. [6ja) for large hx values. One 
can understand this diagram by starting from the dia- 
gram shown in Fig. jsja) and rotating h from the z to the 
x-direction. During this rotation, we expect the phase 
diagram to open at large 1^'s and the CI phases to sep- 
arate and to continuously reduce size. Since there is an 
abrupt simultaneous change in the quantized C± along 
path 2, we have to conclude that all that is left when the 
rotation is completed is a doubly-degenerate line of ex- 
tended states, which terminates abruptly at a critical Wc 
(which is possible if the line is doubly degenerate). For 
W^s lower than Wc^ the mobility spin-gap stays open and 
the change in the quantized C± when moving along path 
2 must be due to the closing of the mobility energy-gap. 
From Figs. |4] and [5| we know that the mobility energy- 
gap stays open when moving along path 1, so the change 
in the quantized C± seen in Fig. [6jb) must be due to 
the closing of the mobility spin-gap. Therefore, as illus- 
trated in Fig. [6]^ a), there must exist a line in the phase 
diagram where the mobility spin-gap closes, and this line 
together with the line of extended states completely en- 
circle the region where C± = ±1. Our final note here is 
that the change in the quantized C± brings measurable 
physical effects even beyond Wc^ where we can still see 
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energy regions with very large localization lengths, levi- 
tating and annihilating each other, in stark contrast with 
what happens in trivial insulators [161 Ull ^ where the ex- 
tended spectrum simply fades away when the disorder is 
increased. 

Summary. In conclusion, we used the spin-Chern num- 
bers to predict the existence of robust bulk extended 
states in TIs with broken TRS and to map out the phase 
diagrams of the models. The predictions were well con- 
firmed by the numerical computations of level statistics 
and localization lengths. Our study settles a long de- 
bated issue, namely, what protects the bulk extended 
states in TIs in the strong disorder regime? The answer 
is: topology alone. In contradistinction with the present 
common belief that the breaking of TRS in QSH insula- 
tors will lead to the sudden and complete localization of 
the bulk states, we found that QSH insulators with large 
Zeeman fields continue to display robust bulk extended 
states. We were able to derive the phase diagrams of the 
systems and, in the process, we discovered that the sys- 
tems can be driven into topological CI phases by disorder 
and Zeeman fields. 
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